NMR Observation of Rattling Phonons in the Pyrochlore Superconductor KOs 2 Oe 
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We report nuclear magnetic resonance studies on the /3- pyrochlore oxide superconductor KOS2O6. 
The nuclear relaxation at the K sites is entirely caused by fluctuations of electric field gradient, which 
we ascribe to highly anharmonic low frequency oscillation (rattling) of K ions. A phenomenological 
analysis shows a crossover from overdamped to underdamped behavior of the rattling phonons with 
decreasing temperature and its sudden sharpening below the superconducting transition temperature 
T c . Suppression of the Hebel-Slichter peak in the relaxation rate at the O sites below T c also indicates 
strong electron-phonon coupling. 

PACS numbers: 76.60.-k, 74.70.-b, 74.25.Kc 



The family of new superconductors ^40s206 (A = K, 
Rb, and Cs) with the pyrochlore structure exhibits many 
exotic properties. For instance, the superconducting 
transition temperature T c increases substantially with 
decreasing radius of A ions (T c = 3.3, 6.3, and 9.6 K for 
A = Cs, Rb, and K, respectively [E & [1, 0] ) and shows 
non-monotonic pressure dependence [5(, in contradiction 
to what would be expected from the change of density of 
states alone. Among this family KOS2O6 is most anoma- 
lous. The electrical resistivity shows strong concave In- 
dependence down to T c [l|, Q , in contrast to the normal 
T 2 -behavior in Rb and Cs compounds @,Sli|. The T- 
linear coefficient of the specific heat 7 = 70 mJ/K 2 mol 
0, 0] is highly enhanced over the value obtained from 
band calculations 10-11 mJ/K 2 mol 0, E3]. 

These anomalies imply large density of low energy ex- 
citations, which was originally ascribed to magnetic frus- 
tration inherent to the pyrochlore lattice, a network of 
corner-shared tetrahedra formed by Os atoms. However, 
the Pauli-like susceptibility Q excludes local moments 
on the Os sites. How the frustrated lattice geometry 
affects properties of itinerant electrons is still an open is- 
sue. On the other hand, there is now a growing body of 
experimental evidence for anharmonic motion of isolated 
alkaline ions in an oversized cage of Os-0 network. The 
specific heat data show existence of low frequency Ein- 
stein modes in addition to the electronic and the Debye 
contributions for all members of AOS2O6 @, HEl]- The 
X-ray results revealed huge atomic displacement param- 
eter for the K sites in KOS2O6 [12J. Recently a second 
phase transition at T p = 7.5 K was observed for high 
quality single crystals of KOS2O6 @- Since the transi- 
tion is of first order and persists into the normal state 
in high magnetic fields, it is likely a structural transition 
associated with the K ion dynamics. Calculation of ef- 
fective ionic Hamiltonian also indicates instability of K 
ions in a highly anharmonic potential [l3j], which might 
be the origin of various anomalies. Oscillations of iso- 
lated ions in a large space have been discussed also in 
other materials such as clathrate [3] and skutterudite 
[IH compounds and commonly called "rattling" . 



In this Letter, we report results of nuclear magnetic 
resonance (NMR) experiments on the 39 > 41 K nuclei (nu- 
clear spin I = 3/2) and 17 nuclei {I = 5/2) in KOs 2 6 . 
We found that the nuclear relaxation at the K sites is 
caused entirely by phonons via the electric quadrupole 
interaction. From a phenomenological analysis we pro- 
pose that rattling phonons undergo a crossover from un- 
derdamped to overdamped behavior with increasing tem- 
perature. Strong electron-phonon coupling is indicated 
by rapid reduction of the damping below T c . On the 
other hand, the relaxation rate at the O sites is sensitive 
to electronic excitations. Substantial damping of quasi- 
particles due to rattling phonons near T c is indicated by 
suppression of the Hebel-Slichter peak at the O sites. A 
part of the 39 K NMR results has been discussed in terms 
of spin fluctuations in a previous publication which 
must now be discarded. 

Two powder samples of KOS2O6 were used in the 
present experiment. The sample A was prepared by heat- 
ing a mixture of KO2 and OSO2 with Ag2 as an oxidiz- 
ing agency The sample B was enriched with 1 O by 
heating KO2 and Os metal in two steps in oxygen atmo- 
sphere containing 45 % 17 0. Both samples show identical 
T c of 9.6 K. Note that the second transition at 7.5 K has 
never been observed in powder samples. The nuclear 
spin-lattice relaxation rate 1/T% at the 39 - 41 K sites was 
measured by the saturation recovery method. Although 
the K sites have tetrahedral (T^) symmetry, hence should 
not suffer quadrupole broadening, the 39 K NMR spectra 
consist of a narrow central line (0.3 kHz HWHM at 8.5T) 
and a slightly broad satellite line (10 kHz HWHM) proba- 
bly due to imperfections. The line shape shows no appre- 
ciable temperature (T) dependence above T c . The entire 
spectrum was easily saturated by rf comb-pulses, result- 
ing in single exponential recovery of the spin echo inten- 
sity. The 17 NMR spectra show a quadrupole broad- 
ened powder pattern. To determine the recovery 
of the spin-echo intensity of the central line as a func- 
tion of time t after the inversion pulse was fit to the form 
[13 M eq - M {UeM-t/Ti) + V exp(-6t/T 1 ) + (1 - U - 
V) exp(— 15i/Ti)} with U and V fixed to the same values 
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at all temperatures. 

The T-dependences of l/(T t T) at the 39 K sites is 
shown in Fig. 1. The inset of Fig. 1(a) displays the 
isotopic ratio 39 Ti/ 41 Ti for the sample A. Surprisingly, 
the isotopic ratio coincides with the squared ratio of the 
nuclear quadrupole moments ( 41 Q/ 39 Q) 2 — 1.48 rather 
than the nuclear magnetic moments ( 41 7/ 39 7) 2 = 0.30 in 
a wide T-range 6-100 K. Thus the relaxation is entirely 
caused by fluctuations of electric field gradient (EFG) as 
opposed to magnetic fluctuations. Since any active elec- 
tronic states should be dominantly s-like at the K sites 
with negligibly small quadrupole coupling to K nuclei, it 
must be phonons that causes the nuclear relaxation. 

The prominent features of the data in Fig. 1 are: (1) a 
peak near 13 K (16 K) for the sample B (A), (2) approx- 
imate constant behavior at high temperatures, (3) rapid 
decrease at low temperatures, (4) a clear kink at the su- 
perconducting T c and sudden decrease below T c . The 
last point indicates strong influence of superconductiv- 
ity on the phonon dynamics, a direct evidence for strong 
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FIG. 1: T-dependence of l/(TiT) at the 39 K sites for the two 
samples at different magnetic fields. The data for a wide T- 
range are shown in (a). The solid line shows the calculated 
result of Eq. (5) as described in the text. The inset shows the 
isotopic ratio of 1/Tj.. The dotted (dashed) line indicates the 
squared ratio of the nuclear quadrupole (magnetic) moments. 
The data in the low T region are expanded in (b). The arrows 
show the superconducting T c at the respective fields. 



electron-phonon coupling. The two samples show slightly 
different results near the peak temperature. Although 
the reason is not understood clearly, this may arise from 
different degrees of imperfection or hydration 0] . 

For nuclei with spin 3/2, the transition probability W q 
between two nuclear levels \I Z = m) and \I Z = m ± q) 
(q = 1 or 2) due to quadrupole coupling is given by the 
correlation function of EFG as [18f 



W q = 



_L/eQ\ 
12 V h J 



([V +q {t),V- q {Q)]W q " Lt dt, (1) 



where ojl is the nuclear Larmor frequency and [A, B] 
stands for (AB + BA)/2. The irreducible components 
of the EFG tensor are defined as V±\ = V xz ± iV yz , 
V± 2 = {V xx - V vy ± 2iV xy )/2, where V xy = d 2 V/dxdy 
etc. is the second derivative of the electrostatic poten- 
tial at the nucleus. For simplicity we assume spherically 
symmetric EFG fluctuations. Then W\ = W% = W and 
the relaxation rate is given as 1/Tj = 2W. 

In diamagnetic insulators phonons are usually the 
dominant source of relaxation for quadrupolar nuclei. 
The ionic motion modulates the EFG, 



y± q = y±q,o 



9V± lu , l d2v ±i u 2 
du 2 du 2 



(2) 



where u is the ionic displacement from the equilibrium 
position expressed as a linear combination of the phonon 
creation (a^) and annihilation (a) operators. The second 
term inserted into Eq. (1) leads to the direct process [ljj ]. 
which is expressed by the one phonon correlation function 
or the imaginary part of the susceptibility, 



W g oc k B T 



Imx(gcJ L ) 
qfuvL 



(3) 



provided qhwL <S k B T . For harmonic phonons with infi- 
nite life time, the contribution from the direct process is 
negligible since u>jj (~ 10 MHz) is many orders of magni- 
tude smaller than the typical phonon frequency and the 
phonon density of states at u>l is practically zero. 

On the other hand, the third term in Eq. (2) leads 
to the two phonon Raman process, in which a nuclear 
transition occurs by absorbing one phonon and emitting 
another phonon [18j. Usually this is by far the domi- 
nant process expressed as 1/Ti oc J{p(uj)} 2 n(uj)(n(uj) + 
l)\A(uj)\ 2 duj, where n(cu) is the Bose factor, p(oj) is the 
phonon density of states, and |^4(w)| is the transition 
matrix element. One can see 1/Xi oc T 2 at high temper- 
atures where n(u)) ~ T/uj. At low temperatures, 1/Ti 
decreases monotonically and approaches either to T 7 - 
dependence for acoustic phonons or activated behavior 
for optical phonons [l^ . Our results in Fig. 1 are, how- 
ever, in strong contradiction with such behavior. 

Since dominance of phonons for nuclear relaxation is 
extremely rare in metals, we suppose that rattling makes 
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a special case for KOS2O6. The distinct features of rat- 
tling are its low energy scale and strong anharmonicity, 
which should lead to damping. We expect that damping- 
induced broadening of the phonon spectrum peaked at a 
very low frequency, will bring substantial spectral weight 
at zero frequency so that the direct process (Eq. 3) be- 
comes important. The peak in \/(T\T) is then natu- 
rally explained as follows. Let us assume a low frequency 
Einstein mode at u>q. At low temperatures, the spec- 
trum Im%(o;)/£j is sharp, therefore, the weight at ojl ~ 
should be small. As the spectral width increases with in- 
creasing temperature, the weight at ui = will grow and 
take a maximum when the width becomes comparable to 
loq. Further broadening, however, will spread the spec- 
trum over higher frequencies and reduce the weight at 
lu = 0. (See the inset of Fig. 2.) Thus the peak in 1/{T X T) 
can be reproduced if the rattling undergoes a crossover 
from underdamped to overdamped behavior with increas- 
ing temperature. The contribution from the Raman pro- 
cess should be also influenced by damping [19(. However, 
we could not reproduce the peak in 1/(T\T) by the Ra- 
man process without fine tuning the T-dependence of the 
damping. We consider that the Raman process is unim- 
portant. The details will be discussed elsewhere [iol |. 

In order to quantify the above picture, we employ 
a suitable phenomenological model describing strongly 
damped oscillators. A commonly used expression is 
the damped harmonic oscillator model, Imx(w)/w = 

t/{(luI~lu 2 ) 2 +t 2 lu 2 }, which lcads t0 i/(r lT ) by 

taking the limit lul — > in Eq. (3). However, this model 
is unable to reproduce the observed peak for any mono- 
tonic change of the damping T. In the damped harmonic 
oscillator model, only the friction or the momentum re- 
laxation is considered as the sources of damping. It was 
argued, however, that in the collision dominated regime 



one must also consider the jump of the spatial coordinate 



One then obtains 
which should be 



or the amplitude relaxation [21 




the van Vleck-Weisskopf formula [2 
suitable for strongly damped anharmonic phonons, 



Imx(w) 



17 w 



(lu + lu ) 2 +T 2 {u-u a ) 2 +T 2 ' 



(4) 



From Eqs. (4) and (3) we obtain in the limit ojl — * 0, 



1 



1 



(5) 



where 7 = T/luq. For a fixed value of loq, this formula 
has a maximum at 7 = 1, which corresponds to the ob- 
served peak. By identifying the T-ranges below or above 
the peak temperature with the underdamped (7 < 1) 
or the overdamped (7 > 1) regions and assuming that 
luq is constant, we extracted the T-dependence of 7 from 
the experimental data of l/(TiT) as shown in Fig. 2. 
The result above T c is well represented by an activation 
law 7 = 7 exp(— E/ksT) with 70 = 5 and E = 20 K 
(the solid line). Thus our analysis predicts extremely 
overdamped behavior at high temperatures. Note that 
the phonon frequency u>q cannot be determined from our 
analysis. The solid line in Fig. 1 shows the T-dependence 
of l/(TiT) calculated from Eq. (5) and this activation 
law. A possible interpretation of the activated behav- 
ior is that r represents the hopping frequency of the K 
ions among different metastable positions in highly an- 
harmonic potential [13] and the activation energy corre- 
sponds to the potential barrier between them. 

We have also measured the relaxation rate for 85 Rb and 
87 Rb nuclei in RbOs2C>6 and separated the phononic and 
magnetic contributions poj . The phonon contribution to 
l/(TiT) in RbOs 2 06 does not show a peak, indicating 
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FIG. 2: The normalized damping 7 = F/uuo is plotted against 
1/T. The arrows correspond to T c at 5 T and 8.5 T. The 
inset shows the phonon line shape Imx(ai)/oi calculated from 
Eq. (4) for several values of 7. 
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FIG. 3: T-dependence of l/(TiT) at the 17 sites at 2 T (solid 
circles) and 8.5 T (open squares). The arrows indicate T c at 
the respective fields. The inset shows the fitting to a power- 
law 1/T oc T a with a = 4.7 (solid line) and an activation law 
1/T tx exp(A/T) with A = 17 K (dashed line) for the data 
at 2 T below T c . 
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underdamped behavior in the whole T-range below the 
room temperature. The T-dependence of 7 in RbOs206 
obtained from similar analysis is also compatible with an 
activation law with E=12 - 20 K, although the accuracy 
is limited. 

A remarkable feature of the results in Fig. 2 is the sud- 
den reduction of 7 below the superconducting T c . This 
indicates that the damping of rattling phonons near T c 
is primarily caused by electron-phonon interaction and 
the phonon life time is enhanced by opening of the su- 
perconducting gap. Therefore, luq should be smaller than 
2 A ~ 50 K [7j, where A is the superconducting gap. 

We now discuss the results at the O sites. Figure 3 
shows the T-dependence of 1/(T]T) at the O sites at 
the fields of 2 T and 8.5 T. Unlike the result for the K 
sites, 1/(T]T) is only weakly T-dependent above T c , ap- 
proximately consistent with the Korringa relation. Thus 
the spin dynamics of conduction electrons are probed at 
the O sites. This is consistent with the Os-5c?(t2 a ) and 
0-2p hybridized nature of the conduction band [H], [l3| . 
However, the weak T-dependence of l/(TiT) in the nor- 
mal state is still anomalous and suggests development of 
modest spin correlation. 

Below T c , it is known that 1/(T]T) shows a peak 
near T = 0.9T C (the Hebel-Slichter peak [H) in weak 
coupling superconductors with an isotropic gap. This 
is not the case for the data in Fig. 3; 1/(T\T) de- 
creases gradually below T c , indicating that the Hebel- 
Slichter peak is strongly suppressed. At lower tempera- 
tures, T-dependence of can be fit either by a power 
law 1/Ti oc T Q with a ~ 5 or by an activation law 
1/Ti oc exp(-A/T) with A = 17K. Since the upper 
critical field H c2 is about 30 T at T = [H, , pair 
breaking effects of magnetic field should not be strong at 
2 T [27|. 

The behavior at low temperatures rules out highly 
anisotropic gap structure such as line nodes, which leads 
to T 3 dependence of \jT\. This is consistent with 
the field independent thermal conductivity showing an 
isotropic gap [28[. Therefore, suppression of the Hebel- 
Slichter peak cannot be ascribed to a highly anisotropic 
gap. It appears more likely due to inelastic scattering of 
quasiparticles by phonons, which broadens the supercon- 
ducting density of states as was proposed for the case 
of TlMoeSe7.s [29]. This has been also predicted by 
numerical calculations based on the Eliashberg theory 
30l 31 1 . The calculation by Akis et al. 31], for exam- 
ple, shows that the Hebel-Slichter peak disappears when 
kBT c /huJi n is larger than 0.2 ~ 0.3, where u>i n is the 
Allen-Dynes parameter representing the typical phonon 
frequency. This is compatible with our previous conclu- 
sion that the frequency of the rattling phonon is smaller 
than 2A ~ 5ksT c . Whether the rattling works as a pair- 
ing mechanism or acts as a pair breaker is still an open 
question. 

In conclusion, our analysis of the l/(TiT) data at the 



K sites indicates a crossover of rattling from the low T 
underdamped to high T overdamped behavior near 13 
- 16 K. The rapid decrease of l/(TiT) at the K sites 
below T c indicates strong electron-phonon coupling and 
enhancement of the lifetime of rattling phonons. Hence 
its frequency must be smaller than 2A. The results at the 
O sites rule out strongly anisotropic gap. The absence of 
the Hebel-Slichter peak at the O sites is also compatible 
with strong electron-phonon coupling. 
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